Adsorption of CO2 and CH 4 and their mixtures in gas hydrates. 
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We report results from grand-canonical Monte Carlo simulations of methane and carbon dioxide 
adsorption in structure si gas hydrates. Simulations of pure component systems show that all 
methane sites are equivalent, while carbon dioxide distinguishes between two types of sites, large or 
small. The adsorbed mixture can be regarded as ideal, as long as only large sites are occupied. A 
strong preference is demonstrated for methane, when the smaller sites become filled. 

The molar heat of adsorption of methane decreases with composition, while the molar heat of 
5— l ' adsorption for carbon dioxide passes an extremum, essentially in accordance with the observation on 

the site sizes. The Helmholtz energies of the hydrate with CO2 -CH4 gas mixture for temperatures 
■ between 278 and 328 K and pressures between 10 4 and 10 9 Pa indicate that certain mixtures are 

more stable than others. The results indicate that a thermodynamic path exists for conversion of a 
pure methane hydrate into a pure carbon dioxide hydrate without destroying the hydrate structure. 
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I. INTRODUCTION 



Understanding the possibility to exchange CH4 with CO2 in hydrates is of great importance to the oil industry, 
and has been the focus of several studies |27H3l| . By studying the radial distribution function of a fully occupied si 
hydrate using molecular dynamic simulations, Geng et al |32| have recently suggested that a CO2 + CH4 mixture in 
I the hydrate lattice can be more stable than a hydrate with pure components, CH4 or CO2. The present work aims to 
find conditions for which this is true, using Grand- Canonical Monte-Carlo simulations (GCMC) in the study of the 
adsorption of pure CO2 and CH4, as well as mixtures of CO2 and CH4 in an si hydrate. 

Monte-Carlo simulations of this type are excellently suited for this purpose, because they allow one to study 
also metastable configurations, in particular, non- fully occupied hydrates [33l43a |. Molecular simulations allows one 
to study the conditions of hydrate formation which are complicated to work with experimentally. For a particular 
• configuration of a hydrate with adsorbed gas(es) one can use the adsorption isotherms to obtain the Helmholtz energy, 
l—— ') which can be used to give information about relative stabilities. GCMC simulations have been used extensively to 
study adsorption processes [HI, H3| and we shall take advantage of the method development described earlier [35H3"" 1 
We shall use the chlatrate si structure, which forms at typical reservoir conditions (for example, pressures up to 10 
\Q . bar or 10 9 Pa, which corresponds to reservoir depths up to 10 km [Hj]). The structure, which is well established 39], 
t^- ' has 8 cages per unit cell. Interaction potentials for water are available through a series of systematic investigations 
[ioj . We will use a rigid water model and the rigid framework model for all of the calculations. The typical behavior 
of the adsorption isotherms remain essentially unaltered when the hydrate framework model becomes mobile (and 
, therefore is not presented here). 

The Helmholtz energy difference is negative for a spontaneous process at constant volume and temperature. The 
I process of practical interest in this study can be written schematically as the reaction 
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N C o 2 + N CHi + N H . 2 o ^ (N C o 2 + N GHi ) ■ N H . 2 o (LI) 

The Helmholtz energy difference per mole of formed unit cell in this reaction is denoted by A r F. The long-term aim 
of the work is to obtain the values for A r F which would decide on whether the hydrate formation is favorable or not. 
This paper concerns a first step in this direction. We will study the variation in the Helmholtz energy of the hydrate 
with a mixed gas adsorbed, (Nco 2 + ^ch 4 ) • Nh 2 o, compared to a metastable hydrate with empty cages Nh 2 o- 
This difference is denoted by AF in this paper. Since, adding a gas molecule in a cage stabilizes the hydrate, this 
Helmholtz energy difference will normally decrease with the loading. However, if a gas hydrate with two components 
adsorbed is less stable than a single component one, the variation in the Helmholtz energy would also increase with 
loading. The main purpose of the paper is to show this, and discuss how this process can be understood in terms of 
molecular behavior. The quantities A r F and AF are related, but we return to the relationship and a discussion of 
A r F in a subsequent paper. 

The paper is organized as follows. In|H]we discuss the basic thermodynamic relations which allow us to calculate 
the Helmholtz energy and the molar heat of adsorption as a function of loading. The details of the simulations are 
specified in IIII1 In IIVI we provide our results. We report adsorption isotherms for the single component hydrates and 
the mixture of CO2 and CH4 in the hydrate, as well as the selectivity data for the mixture adsorption. Furthermore, 
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we calculate the Helmholtz energy difference AF of the above hydrates as a function of loading. Finally, we compare 
the data for the molar heat of adsorption with the structure of the hydrates. An overall discussion and conclusion are 
given in [V] 



II. THERMODYNAMICS OF ADSORPTION 



A. Helmholtz energy 



We consider the framework of water molecules in the si clathrate structure. In simulations we can control whether 
this framework is flexible or rigid. In both cases we specify and keep fixed the lattice parameters of the framework, 
so the volume of the system V is constant. We also keep the temperature of the system T fixed. 

The Helmholtz energy F of the clathrate with given numbers of adsorbed guest molecules (a loading) is given by 

n 

F = —pV + (j, w N w (11.1) 
i=i 

Here p is the pressure in the system, ^ is the chemical potential of the i-th guest component and N is the number 
of adsorbed molecules of component i, while fj, w and N w are the corresponding quantities for water. The number of 
compounds is n. For pure components, n = 1, for a mixture of methane and carbon dioxide n — 2. The number 
of water molecules is fixed, while the number of guest molecules vary. The chemical potentials, and therefore the 
Helmholtz energy, depend on the loading. 

We would like to calculate AF(N^y) = F(N^) — F(0), the difference in the Helmholtz energy of a filled clathrate 
F(N {l} ) = F(N w ,Ni,--- ,N n ) and the Helmholtz energy of the empty clathrate F(0) = F(N w ,0,--- ,0) = -pV + 
N w . We refer to this difference as the Helmholtz energy of the mixture of gases in the hydrate. It follows from III. II 
that 

n 

AF(N {i} ) = - [ P (N {i} ) - p(0)] V + M^«) ~ MO)] N w +Y, IH(N {iy ) N (11.2) 

i=l 

where p(N^y) is the hydrostatic pressure at the given volume and loading, while p(0) is the hydrostatic pressure at 
the same volume and zero loading. According to the Gibbs-Duhem relation we have 

n 

SdT-Vdp + N w dn w +^2N l dfi l =0 (11.3) 

i=l 

Integrating the equation at a constant temperature, volume, and a constant number of water molecules from zero 
loading to a given loading we obtain 

- V [p(N {i} ) -p(0)] +N W [fx w (N {i} ) - MO)] / NidiH(N {i} ) = (II.4) 







which, after substitution into III. 2 1 and integration by parts, gives [35| 

n n pnU N w) 

AF(N {1} ) = Vi(N {i} )Ni - / N * d ^( N H}) 

i=i i=i Jo 

where fx* is the specified chemical potential of each of the component. The right hand side of this expression can be 
easily obtained from GCMC simulations, in which one specifies the chemical potential of the adsorbing components 
and measures the number of adsorbed molecules. Furthermore, this form is convenient for integration, because at low 
loadings the slope of the function Ni(/J,i) approaches zero. The integration in III.5l over the chemical potential of one 
component is done by keeping the chemical potentials of all other components fixed. For a two-component system 
this means that integration is performed first along the path of \ii = —oo from fi\ = —oo to fi\ = fxl and then along 
the path fii = /i* from (j,2 = — oo to ^ = A*2- 

The chemical potential of the i-th component directly follows from its fugacity ft [U [42| : 

IH = -k B T\n(R qi ) + k B Tln(fi) (II.6) 
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where the first term in lll.6l is the contribution from the ideal gas. Here qi = Z- lnt ^l '(iV^A?) where Z lnt ,i is the partition 
function for internal degrees of freedom. If the molecules are rigid (which is the case for our simulations, see 
^int,i has contributions only from rotation of the molecule. At high temperature Z lnt ;\ = T/T r , where T r is the 
characteristic temperature of the rotational degrees of freedom. Furthermore, A< = HNa/ ^/ZirMiRT is the thermal 
de Brogile wavelength of component i, Mi is the molar mass, and h, Na, R, ks are Plank's constant, Avogadro's 
number, the universal gas constant, and Boltzman's constant, respectively. 



B. The heat of adsorption 



The heat of adsorption of gas(es) into the hydrate at constant volume is equal to the internal energy change of the 
reaction at constant volume. From the internal energy of reaction (jl.ip . we define the partial molar energy for 
adsorption of a component i; 

ArUl . . ( a JE^E!l) (n .r) 



\ dNi J Nj,V,T V 9Ni / N j: V,T 

where U h is the internal energy of the hydrate phase and U 9 is the corresponding value for the gas phase. For a guest 
molecule in the hydrate, both energies have constant contributions from the translational energy, ^ksTN^ where 
Si is the number of degrees of freedom of the component, and these contributions cancel in the difference 43]. The 
remaining contribution is due to the difference between the configurational parts of the internal energy U£ and U§. 
In the literature it is common to consider the gas to be ideal, so that it acts as a reference state to an adsorbed phase 
(3(||43j]. Therefore, Uf c — NiRT. This defines the partial molar heat of adsorption q as 



- RT (II. 

N U T,V 



In this study we are interested only in the hydrate phase. The quantity —qi defined by 111.81 is then a property 
of a hydrate phase with reference to an ideal gas phase. This is typical in the simulation community and conforms 
with the discussion in III Al to describe the hydrate phase only [36j . We shall use 111.81 for hydrates filled with single 
components, and calculate the molar heat of adsorption of CH4 and CO2. 

The derivative of the internal energy of the hydrate with respect to the number of adsorbed molecules can be 
calculated in a number of ways [36j. In GCMC simulations, one calculates the molar heat of adsorption (111.81) from 
the fluctuations of the thermodynamic quantities as [13, IH| 

(5U?6N) 

q - RT -lsNSNy (IL9) 

where (SX SY) = (XY) — (X)(Y), and () denotes an ensemble average in the grand-canonical ensemble. 



III. SIMULATION DETAILS 



We perform grand-canonical Monte Carlo (/jXT) simulations of CH4 and CO2 in a 2x2x2 unit cell of si hydrate, 
which has 64 cages, with a lattice parameter 12.03 A [39| . The number of running cycles was varying from 500 
to 500000 in order to achieve the desired accuracy. The system became equilibrated rather quickly, typically after 
300-500 cycles. The number of MC moves per cycle is equal to the number of particles of each component in the 
system, with a minimum of 20. 

In GCMC simulations one specifies the chemical potential of a component and calculates the average number of 
particles which correspond to this chemical potential. For this a number of Monte Carlo moves is performed, typically 
attempts to displace a particle from one position in the box to another and to exchange particles with the reservoir. 
In the semi-grand canonical MC simulations, we specify N w and /, or fii, i = l..n and we find the loading of the guest 
components Ni and accordingly the Hclmholtz energy of the gas mixture in the hydrate. In our GCMC simulations 
the molecules were allowed to change their position and orientation. They were also subjected to Regrow, Swap, 
and Identity change [44[ MC moves. Furthermore, a series of NVT simulations has been performed to analyze the 
distribution of molecules between the cages. In this case the number of adsorbed molecules was specified, while the 
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chemical potential was adjusted during the simulation. The guest molecules were allowed to change the position and 
orientation as well as regrow. 

The oxygen positions of water molecules were taken from crystallographic data [39| . The orientation of hydrogen 
atoms was chosen random in accordance with the Bernal- Fowler rule [451 ] . Initially all the hydrogen atoms were 
assigned randomly to each of 4 sites of the water molecule (2 sites from the hydrogens and 2 sites for the hydrogen 
bonds to the neighboring molecule). Then a short Monte Carlo procedure has been performed to displace the hydrogen 
atoms between the sites randomly, in order to satisfy 2 conditions: i) each oxygen atom has only 2 hydrogen atoms, 
and ii) there is only 1 hydrogen atom between each 2 oxygen atoms. The third requirement, that the number of the 
hydrogen atoms of the corresponding Wyckoff type is given, was satisfied automatically after the above procedure. 

The water model is TIP5PEw [46(. We performed simulations with both immobile and mobile water molecules, so 
in the latter case they were allowed to change their position and orientation. During the simulations we have found 
that behavior of the system with mobile and immobile water molecules is similar, while the fluctuations were larger for 
the case of mobile molecules. In the analysis of the results we therefore present the data for a rigid water framework. 
The description of the guest molecules was taken from the TraPPE force field [47| . 

We performed simulations for the range of the temperature from 278 K to 328 K. The adsorption isotherms were 
qualitatively the same for all the temperatures in this range. We therefore focused on the temperature of 278 K. 

The adsorbed components were assumed to be in equilibrium with an imaginary gas reservoir. Specifying the 
chemical potentials of the components in this gas determines therefore the chemical potentials of the components, 
adsorbed into the hydrate. This procedure is typical in the molecular simulations of adsorption [48(. The chemical 
potentials of the components can be specified in two ways. First, the pressure of the gas in the reservoir and its 
composition can be specified. The Peng-Robinson equation of state is used to calculate the chemical potential or 
the fugacity of each component. This method was used to obtain single component adsorption isotherms and x-y 
diagrams for mixtures. The pressure varied in the range between 10 4 Pa and 10 9 Pa in these simulations. Second, 
one can specify the fugacities of the components directly using the Lewis- Randall rule [42| . According to the Lewis- 
Randall rule, the fugacity of each component equals the product of its mole fraction and the pure component fugacity 
at the same temperature and pressure of the mixture. This method was used to obtain the adsorption isotherms 
for mixtures. This allowed us to perform the integration of 111.51 along each chemical potential, keeping the chemical 
potentials of the other components constant. The fugacity of each of the components varied in the range between 10 4 
Pa and 10 9 Pa. 



IV. RESULTS OF THE SIMULATIONS 



A. Single component adsorption isotherms 

The single component adsorption isotherms for CO2 and CH4 at 278 K are shown in [TJ We see how the hydrate 
loading, N, varies with applied pressure p on a si unit cell with 8 cages. The pressure axis is given in a logarithmic 
scale. The error bar is less than the size of the symbol, for a hydrate with immobile water molecules. The adsorption 
isotherms of CO2 and CH4 for the case of mobile water molecules are also similar (not shown). 



At first sight, the adsorption isotherms in 1(a) may look like Langmuir adsorption isotherms, however the fitting of 
the one- or two- site profiles to such isotherms is not satisfactory. Still, the CO2 isotherm reveals a two-site adsorption 
behavior. This can be explained by the structure of si clathrate. A unit cell of an si hydrate has 6 large and 2 small 
cages. Being a relatively large molecule, CO2 tends to occupy large cages first. Small cages start to be filled only after 
the large cages have been occupied. This can be also seen on[2j This figure shows the distribution of CO2 molecules 
between several cages, obtained from a NVT simulation. The positions of these molecules from 500 snapshots were 
combined and drawn on a single figure. When the loading of CO2 is small (N = 11), the loading of the large cages 
is much larger then the population of small cages. At intermediate loading (N = 52), which corresponds to the 
plateau in the CO2 adsorption isotherm, some small cages start to be filled, while the other small cages remain almost 
unoccupied. Finally, at high loading (N = 66), all the cages are occupied uniformly. 

The cell with 2x2x2 units of si clathrate has 64 cages. Having 66 molecules adsorbed into clathrate means that some 
cages contain 2 CO2 molecules. This happen at large pressures only. Double occupancy of a cage can be explained 
by the linear form of CO2 molecule, which allows special orientation in large cages. 

The CH4 isotherm, on the other hand, reveals a one-site adsorption behavior. The reason for this is that CH4 
molecules are rather small and have almost no preference between large and small cages. However, the cage size is 
not completely irrelevant for the CH4 distribution: otherwise the adsorption isotherm would be a perfect single-site 
Langmuir curve. The typical distribution of CH4 molecules is shown on[3J The snapshots were obtained in the same 
way as those on[5]from fjYT simulations at different gas pressures. This figure shows that CH4 fills the cages uniformly 
showing almost no preference to large or small cages, in contrast to CO2. 
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Furthermore, we computed the adsorption isotherms for different temperatures in the range between 273 K and 328 
K. The results of simulations are presented on [U One can see, that increasing the temperature leads to the decrease 
in the amount of adsorbed molecules both of CO2 and CH4. Nevertheless, the hydrate becomes completely filled 
at pressures about 10 9 Pa. It is interesting to analyze the distribution of the CO2 molecules in cages with respect 
to the temperature. Recall that if the loading is less then 6, then mostly the large cages are being filled, while the 
higher loading correspond to filling also small cages. We see from 4(a) that increasing the temperature makes the 
isotherm to cross the line N = 6 at higher pressure. Furthermore, the plateau at N = 6 becomes less pronounced 
when the temperature is increased. This shifts the type of adsorption isotherm from 2-site to 1-site, which means 
that the difference between cages becomes less, when the temperature is increased. The CH4 adsorption isotherm is 
qualitatively the same at different temperatures. Increasing the temperature only increases the pressure at which the 
desired loading is achieved. 



B. Mixture adsorption isotherms 

Isotherms for hydrates loaded with mixtures of CO2 and CH4 are shown in [5l This figure shows the CO2 mole 
fraction in the hydrate as a function of the gas pressure and the corresponding CO2 mole fraction in the gas phase. 
One can identify two regions in the figure, with the pressure above approximately 10 7 Pa and below 10 7 Pa. As one 



can see from 1(a) this is approximately the pressure where the adsorption isotherm of CO2 reaches the plateau. In 
both regions the molar content of CO2 in the hydrate increases gradually when the content of CO2 in the gas phase 
increases. When the pressure is below 10 7 Pa the CO2 molecules fill the large cages while the CH4 molecules fill both 
the small and the large cages. Thus, there is no preference between large and small cages. Therefore, the mixture 
composition in the hydrate is approximately the same as in the gas phase over the whole range of compositions, see 
[fJJ At approximately 10 7 Pa the large cages are all filled and only small cages are available. CO2 molecules cannot 
compete with CH4 molecules for these cages, since the former ones are larger than the latter ones. It is mostly CH4 
molecules which occupy the small cages in the mixture, see[(JJ 

[5] compares the results of the mixture loading with ideal adsorption solution theory (IAST) [49| . The concept of 
ideal solution for the adsorbed components is analogous to the Raoult's law for vapor-liquid equilibrium. It holds well 
for mixtures of similar components which are adsorbed on similar sites. For segregated systems (when one component 
adsorbs on one type of sites while the other component adsorbs on the other adsorption sites) it is known that it does 
not work [5(| [Hj]. As one can see from[f51 the mixture can be considered ideal at small pressures, which are below the 
pressure when CO2 loading reaches the plateau. At these pressures CO2 and CH4 molecules occupy the large cages 
equally likely. Even though CH4 and CO2 molecules are different, the large cages are so large, that they show no 
preference to the components they would adsorb. At pressures larger than 10 7 Pa we observe preferential adsorption 
for CH4 molecules, so the mixture cannot be considered ideal. At these pressures the small cages are being filled, 
so it is natural that they prefer small CH4 molecules rater than large CO2 molecules. It becomes clear that at high 
pressures IAST does not hold. 

We next consider the dependence of loading on the fugacity of each of the component, which is shown on [71 We see 
that the loading of each of the component increases proportionally to its fugacity. Furthermore, the increase of the 
fugacity of one of the components complicates the adsorption of the other component and vice versa. This explains 
the dark triangular region in the bottom-right and the top-left corners of the figure for CO2 and CH4 respectively. 
Furthermore, one can see the gradual filling of the hydrate with the increase of the total fugacity Jch 4 + fco 2 which 



is related to the gas pressure. This behavior is consistent with the adsorption isotherms on 1(b) 



C. The Helmholtz energy of the pure gases and gas mixtures in hydrates 

We first report the Helmholtz energy for the single component hydrates. [5] shows this Helmholtz energy for CO2 
and CH4 hydrates. One can see, that the more molecules are adsorbed, the lower is the Helmholtz energy per mole. 
This means that the hydrate is becoming relatively more stable as the number of molecules increases. This is well 
known (28|, because empty cavities are not stable and guest molecules are required to stabilize them. Furthermore, 
one can see that the Helmholtz energy for the CO2 hydrate is lower than that of the CH4 hydrate over the whole 
range of loadings. This indicates that the CO2 hydrate is more stable than the CH4 hydrate. 

Furthermore, we consider the temperature dependence of the Helmholtz energy of a single component hydrate. The 
computed data in the range of temperatures between 273 K and 328 K are presented on[9] One can see, that variation 
of temperature does not affect the Helmholtz energy at all. This means, that the change of the entropy with respect 
to the change in loading is zero, AS — (dAF/dT) v N = 0. This indicates that the adsorbed molecules in different 
cages do not interact. 
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We next consider the Helmholtz energy of the gas mixtures in hydrate. [TOl shows how the Helmholtz energy depends 
on the fugacities of each of the components. We see that there is a certain region where this energy is lowest, so 
the hydrate with gas mixture would be more stable at these composition. At low fugacities there are not that many 
molecules adsorbed, which means that many cages are empty. This makes the gas hydrate structure unstable, which 
explains the dark region on the figure in the bottom-left corner. If both fugacities are high, both CO2 and CH4 
molecules prefer to occupy the cages. Competition for an available cage makes the hydrate structure less stable, 
as well. This is indicated by the dark region on the hgure in the top-right corner. Finally, if only one of the two 
fugacities is large, while the other is low, only the high fugacity molecules tend to fill the cages and the hydrate 
structure becomes relatively more stable. 

These considerations of stability are not conclusive, however. As explained in the introduction a more complete 
set of thermodynamic data is needed to evaluate the possibility of hydrate formation. Namely, the values of the 
Helmholtz energy of different phases which can be formed of CH4, CO2 and water are required. We shall return to 
this point in our next work. 

It is still interesting to plot [10] in a p-y scale, where p is the gas pressure, which is in equilibrium with the guest 
molecules in the hydrate, and y is the mole fraction of CO2 in the gas phase. The distribution of the Helmholtz energy 
of the hydrate with gas mixture for this case is shown in 11(a) The pressure and the mole fraction were calculated 
from the fugacities of the components using the Lewis and Randall rule [42j . We see that the hydrate is more stable 
at high pressures (i.e. at high loadings) if pure gas is being adsorbed. In contrast, if we have a mixture, the highest 
stability is achieved at moderate pressures. 



11(a) shows that there is a region in the pressure-composition diagram, which connects a relatively stable methane 



hydrate with a relatively stable carbon dioxide hydrate. While at low and high pressures the hydrate with gas mixture 
is unstable, it is almost as stable as cither single component hydrate at moderate pressures. This is very promising, 
since one of the goals of our research is to find a path on a phase diagram to convert a pure CH4 hydrate into a pure 
CO2 hydrate without destroying it s struc ture. In a pursuit of this path it will therefore be interesting to focus on the 
lightest region on the diagram. On 11(b) this path is drawn schematically. Note that following this path implies that 
one does not meet any energy barrier: the value of AF along it is almost everywhere constant. In contrast, going 
from pure CH4 hydrate to pure CO2 hydrate at constant pressure around 10 s — 10 9 Pa requires crossing an energy 
barrier of 50-100 kJ per mol of unit cell. 

We see that this path goes through a region where a mixture of CH4 and CO2 is adsorbed into a hydrate. Namely, 
one has first to decrease the pressure of pure CH4 hydrate down to approximately 10 7 Pa. At this pressure there is 
enough space in the large cages for CO2 molecules to go into them. Next, one has to increase the content of CO2 in 
the surrounding fluid. This would lead to the filling of the large cages with CO2 molecules while the CH4 molecules 
will be expelled out of the hydrate. When there is no CH 4 in the mixture, the small cages are empty at the pressure 
10 7 Pa, and this is the time to start filling them in with the CO2 molecules. This is the last part of our path which 
ends up with a pure CO2 hydrate. 



D. Heats of adsorption 

[T2l shows the partial molar heats of adsorption, 111.91 of the single component hydrates obtained from GCMC 
simulations. At high loading (N sa 7) the calculated value of the partial molar heat of adsorption starts to fluctuate. 
This is expected for GCMC simulations, since in the dense phase which is realized at high loadings, the fraction of 
accepted trial MC moves is quite low. This means that that we shall trust the results at small and moderate loadings 
only (not more then 6-7 adsorbed molecules per unit cell). Both sets of data agree with the findings above. 

We see in [T^] that the partial molar heat of adsorption of CH4 decreases almost linearly with the loading over the 
whole range of loadings. We also know that the change of the entropy AS with respect to the change in loading (see 



II. ip is zero. The single-site adsorption isotherm shown on 1(a) for methane can thus be explained by 1121 alone. In 



contrast, for carbon dioxide one can distinguish the two regions for the partial molar heat of adsorption of CO2 ■ 
At small loadings the partial molar heat of adsorption decreases, similar to that for CH4. At this moment the large 
cages of the hydrate are being filled in. At moderate loadings the partial molar heat of adsorption has a tendency 
to increase, however. This corresponds to the filling of the small cages of the hydrate. As the CO2 molecules prefer 
not to be in the small cages, more and more energy is required to get them into hydrate at high pressures. Knowing 
that AS = also here, the shape of the curve in [12] can explain the adsorption isotherm for CO2 shown on l(a)| 



Indeed, the profiles on[T2]at small loadings are not exactly linear, which implies that the adsorption isotherm is not 
Langmuir-like. 

The profile for CO2 in [T^] changes direction not at N = 6 which corresponds to the complete filling of the large 
cages. The minimum is when only 4-5 cages are filled in. This may be explained by a strong interactions of large 
asymmetric CO2 molecule with water cages [3l| . 
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The partial molar heat of adsorption of CO2 is lower than the one of CH4 over the whole range of loadings which 
can be trusted. This supports the preferential filling of the large cages by CO2 molecules. 

Furthermore, we consider the temperature dependence of the partial molar heat of adsorption of a single component 
hydrate. The results of simulations in the range of temperatures between 273 K and 328 K are presented on[5J Within 
the error of simulation one can consider the partial molar heat of adsorption to be temperature independent, like the 
Helmholtz energy, computed in the previous section. 



V. CONCLUSIONS 



In this paper we studied the thermodynamic properties of an si hydrate filled with CO2, CH4 and the mixture 
of those. By performing grand-canonical Monte Carlo simulations, we obtained adsorption isotherms for single 
components and gas mixtures, partial molar heats of adsorption of single components, and the Helmholtz energy of 
a hydrate with the gas mixture referred to an empty hydrate. 

The results show that the state of the hydrate is determined by the relation between the size of the clathrate cages 
and the guest molecules. The si hydrate unit cell has two types of cages, 6 large and 2 small ones. The small methane 
molecules can easily fill both the large and the small cages, while the larger and not spherically symmetric molecules 
of carbon dioxide prefer the large cages. This determines the two regions of the hydrate state with various loading. 

At a small pressure the CO2 molecules tend to fill the large cages only. Every new adsorbed molecule decreases the 
Helmholtz energy of the gas mixture in the hydrate so that the hydrate becomes more stable. This happens until all 
the large cages are occupied. Then the small cages start to be filled. While the hydrate is still becoming more stable 
(the Helmholtz energy is decreasing), the relative increase in stability becomes smaller. 

The methane adsorption reveals almost no preference to cage type, so the partial molar heat of adsorption is always 
decreasing, which means that every new molecule is more welcome than the previous one. The Helmholtz energy is 
always decreasing, which means that the fully occupied hydrate is the most stable one comparing to a partially filled 
hydrate. 

Furthermore, we studied the temperature dependence of the loading, the Helmholtz energy and the partial heat of 
adsorption. While the adsorption isotherm of single-component hydrate shifts when the temperature is changed, the 
dependence of the Helmholtz energy and the partial heat of adsorption on loading is unaltered by the change of the 
temperature. This indicates that the adsorption may be considered as ideal, i.e. the interaction between adsorbed 
molecules in different cages is low. 

The single-component adsorption behavior determines the regime of mixture adsorption. The adsorbed mixture 
may be regarded as ideal, but only when the large cages are being filled. At this point there is almost no difference 
between the behavior of CO2 and CH4 molecules. In contrast, when small cages start to be filled, the hydrate with 
gas mixture reveals strong preference to CH4 adsorption. 

Finally, the diagram obtained for the Helmholtz energy of the gas mixture in hydrate is promising. It may suggest a 
way to convert CH4 hydrates, which are excessively available at the sea bottoms, to CO2 hydrates. This could play a 
role for the question of CO2 storage, without destroying a hydrate. If one considers 11(a) as a starting point, a possible 
path corresponds to going from the top left corner (pure CH4 hydrate, full loading) to the top right corner (pure CO2 
hydrate, full loading). In order to keep hydrate stable, i.e. to not use extra energy to destroy it, one has to go along 



the colored path on 11(b) From the kinetic point of view it may be complicated to perform an exchange of molecules 
in a hydrate, since, unlike zeolites, a hydrate does not have empty channels for molecules to go in. In practice, the 
exchange could be performed in a distorted clathrate structure [52}. Existence of such a path must presently be seen 
as qualitative, as we are lacking information on the thermodynamic properties of the the other phases which can be 
in equilibrium with the hydrate phase, cf. 11.11 Our further efforts are now directed into calculating the phase diagram 
and the remaining thermodynamic data. 
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FIG. 1. Number of adsorbed molecules per unit cell of a si hydrate as a function of the applied pressure (a) and fugacity (b) 
as computed by GCMC simulations at T = 278 K. 
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FIG. 2. Distribution of CO2 molecules in the clathrate cages of a cell with 2x2x2 units of si hydrate at various loading: (a) 
N = 11 molecules, (b) N = 52 molecules, (c) N = 66 molecules. Carbon atoms are colored blue, oxygen atoms are colored red, 
water molecules are not displayed for clarity. 
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FIG. 3. Distribution of CH4 molecules in the clathrate cages of a cell with 2x2x2 units of si hydrate at various loading: (a) 
p = 10 4 Pa, (b) p = 10 5 Pa, (c) p = 10 s Pa. CH 4 molecules are colored blue, water molecules are not displayed for clarity. 



11 



9 \ 

8 
7 
6 

I 5 
S 4 

3 

2 

1 



t ****** 




► 

< - - 
★ 



10' 



V 



10 c 



p/[Pa] 

(a) 



10 c 




***** * 



► 

★ 



• 


273 


K 


• 


278 


K 


• 


283 


K 




288 


K 


• 


293 


K 




298 


K 


♦ 


303 


K 


▲ 


308 


K 




313 


K 


► 


318 


K 




323 


K 


★ 


328 


K 


1 




273 


K 




278 


K 




283 


K 




288 


K 




293 


K 




298 


K 


♦ 


303 


K 


▲ 


308 


K 




313 


K 


► 


318 


K 




323 


K 


★ 


328 


K 



JO 



10 6 10 8 
p/[Pa] 

(b) 



10 



10 



FIG. 4. Number of adsorbed molecules per unit cell of a si hydrate as a function of the applied pressure of (a) CO2 (b) CH4 , 
as computed by GCMC simulations at different temperatures. 
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FIG. 5. Mole fraction of CO2 (color) in the CO2 +CH4 gas mixture in the hydrate at 278 K computed by GCMC simulations 
as a function of the gas pressure and the mol fraction of CO2 in the gas described by the Peng-Robinson equation of state. 




FIG. 6. Mole fraction of CO2 in the gas mixture, y, as a function of its mole fraction in the hydrate, x, at 278 K for 
gas pressures, as computed by GCMC simulations and compared to the predictions of ideal adsorption theory (I AST). 
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FIG. 8. Helmholtz energy of a single component hydrate at 278 K as a function of the hydrate loading per unit cell. 
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FIG. 9. Helmholtz energy of a single component hydrate at different temperatures as a function of the hydrate loading per 
unit cell for (a) CO2 (b) CH 4 . 
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FIG. 10. Helmholtz energy difference (color) as computed by III. 51 of the CO2 +CH4 gas mixture in the hydrate at 278 K as a 
function of the fugacity of each of the component. 




FIG. 11. (a) Helmholtz energy difference (color) as computed by |H.5l of the CO2 +CH4 gas mixture in the hydrate at 278 K as 
a function of the gas pressure and mol fraction of CO2 in the gas phase, (b) The path on the diagram to convert a pure CH4 
hydrate to a pure CO2 hydrate. 
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FIG. 12. Partial molar heat of adsorption of a single component hydrate as computed by GCMC simulations at 278 K as a 
function of the loading of a single-component hydrate. 
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FIG. 13. Partial molar heat of adsorption of a single component hydrate as computed by GCMC simulations at different 
temperatures as a function of the loading of a (a) CO2 (b) CH4 hydrate. 
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